
file <- commandArgs(TRUE)[1]
scenario <- commandArgs(TRUE)[2]
descr <- commandArgs(TRUE)[3]

d <- read.table(file)
e <- qnorm(0.975)*sd(d$V1)/sqrt(length(d$V1)) 
m <- mean(d$V1)

cat(scenario, m, m-e, m+e, descr, "\n")

